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Alphasat wireframe model (deployed). 
(Photo: ESA) 
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Experiment Goals 


* To assess the impact of rain attenuation and scintillation effects on links 
operating in the Q-band 


* To develop a physical model to improve predictions of atmospheric attenuation 
within the desired spectrum 


¢ To conduct long term site diversity experiments at Q-band in coordination with 
Spino d’Adda site (20km) 


NASA Motivation 


¢ NASA’s allocation of 4 GHz of contiguous bandwidth in the Q-band provides an 
opportunity to meet next generation bandwidth requirements 


¢ Preliminary architecture studies of the next generation TDRS system will require 
higher downlink bandwidths than available in the current Ku-band allocation 


Launch of Alphasat on an Ariane 5, July 2013. 
(Photo: ESA / CNES / ARIANESPACE) 
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Beacon Receivers Disdrometer 
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Drop Size Distribution (DSD) data from a disdrometer can be 
used to derive a scaling factor to scale measured attenuation 
data to other frequencies. 


The NASA-POLIMI Alphasat terminal in Milan has the capability 
to directly measure 20 GHz and 40 GHz attenuation as well as 
concurrent DSD data, and thus is uniquely situated to assess and 
characterize these techniques. 


Site of Study 


Milan, Italy 

45.4787° N, 9.2327°E 
Rainfall: 67 mm / month 
Avg. Humidity: 73% 


Milan, Italy 


Installation Date April 2014 


Ground Latitude 45.4787° N 
Station —_ Longitude 9.2327°E 
— Altitude 138 m 
Name Alphasat 
e Nom. Elevation 35° 
al : 
mh Satellite om. Azimuth 158° 
alpha 
Beacon Freqs. 19.701 GHz 
39.402 GHz 
Politecnico di Milano, DEIB Building 5 


(Photo: Google Earth) 
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Beacon Receivers Optical Disdrometer Weather Instrumentation 


Anemometer: 

Young 05178A 
Temperature/Humidity Sensor: 
Young 41382VC 

Pressure Sensor: 

Young BPV3000 

Tipping Bucket: 

Young 52203 


Antenna Gain | 45.6 dBi (Ka / Q-band) 
ling P J | 0.125 sec (8 Hz) 


Beacon Receiver Design 


NA Downconverter Card 


Ka-Band RF IF &ADC 
19.701 GHz 
SS AL 70 MHz to 455 kHz 


ADC 
NI PCI-5124 


fs = 1.111 MHz 
Q-Band RF : Q N=2" 


39.402 GHz 
LO 
19.203 GHz 


Beacon Receiver Specifications 
Jow Sey ee) )) | 3-step down to 455 kHz 


: BPF BPF IF Amp. 
on Y7areceyahv(=1a-J(e) am (0) | 4-step down to 455 kHz 20.2 GHz 70 MHz +25 dB 

504 K (Ka-band) 

| 720 K (Q-band) 
_ 38 dB (Ka-band) ¢ Common ultra-stable 10 MHz ref. oscillator 

40 dB (Q-band) ¢ T xa temperature controlled within +0.1°C 
=| a Mane * Trerpiate temperature controlled within +0.01°C 
2 e 


Bee ee anne T |, temperature controlled within +0.25°C 
See sel 8 Hz / 1 Hz (averaged) 


Disdrometer Specifications 


Particle 


* The disdrometer utilizes an infrared laser diode to generate a 785 56 W/ t 
mm y 0.75 mm 


nm parallel light beam covering a 4560 mm area swath. 


Infrared Light Beam 


* A photodiode on the receive end monitors the reduction in signal 
strength as precipitating particles cross the beam. From this 
measurement, details on the particle size and fall velocity can be 
derived from the magnitude and duration of the signal attenuation, 
respectively. 


228 mm 


* On-board DSP classifies the type, intensity, quantity, and particle 
spectrum of precipitation over 1 minute intervals, as well as 
estimating rain rate from the DSD. 


¢ Spectrum is classified into 22 diameter bins from 0.125 mm to 8 
mm and 20 velocity bins from 0 to 10 m/s with non-uniform bin 
widths. 


Depiction of disdrometer particle detection. 
(Photo: Thies Clima) 


Specific Attenuation from DSD Data 


Ga a : >» 
Specific attenuation is a function of the Specie Stenuaion 
wavelength, forward scattering coefficient, and 22 
drop density distribution [5]. ne 4.343 x 102 as > Re{S(0)}N(D)AD 
7 
The forward scattering coefficient is Geeeincate ti 
dependent on frequency, radius of the rain drop, if a ae aad pad 
and temperature, and is calculated utilizing the Re{S(0)} Mea rere Cneliieiont 
Bie statiemng moder le) N(D) Drop Density Distribution 
AD Drop Size Interval J 
. 
a Drop Density Distribution ) 
n;(D;)x10° 
A COS I rreerremeraaer 
v(Dj)xXSXTXADj; 
Drop density distribution is calculated from . : , 
the binned DSD data using the relation shown a z ith bin : : 
[7]. nD) Number of Drops with Mean Diameter D, 
v(D)) Terminal Velocity in Still Air 
; 2 
All of these parameters are directly available > a paloa 2) [ALG Nt) 
from the recorded DSD data and _ the T pes Time (60 a 
specifications of the disdrometer hardware, L AD; Bin Width of Drop Size Class J 
except for the terminal velocity, which is 
modelled using the analytical fit of measured Terminal Velocity 


data by Gunn and Kinzer [8-10]. 


v(D;) = 9.65 — 10.3e7 9% 
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Instantaneous Frequency Scaling 
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The prediction tracks well throughout the 
event. 
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Heavy Rain Event (<40 mm/hr) 


Attenuation exceeds the Q-band system 
dynamic range (40 dB) from approximately 
00:45 — 01:00 UTC, while the DSD predicts 
attenuation up to 96 dB during this period. 


The figures above plot examples comparing the measured Ka- and Q-band attenuation to the DSD- 
derived Q-band attenuation calculated by applying the instantaneous scaling factor to the Ka-band 


measurement. 
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Rain Event along Path (<10 mm/hr) 


Attenuation up to 11.6 dB is measured during 
the first half of the event, but no rain is 
measured at the site, indicating that the rain is 
along the path. By the end of the event, the 
rain cell appears to have moved overhead the 
receiver, and the DSD prediction tracks the 
measurement. 
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Mie vs. T-Matrix 


Extinction Cross Section (mm2) 
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Change in Specific Attenuation 


Specific Attenuation Difference in Specific Attenuation (T - Mie) 
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Scattering Model Accuracy 


. Measured Attenuation 
Measured Attenuation 


2015-04-27 
2016-04-23 - z : : Measured Attenuation 
; 1 1 jG 2016-05-19 
35H Ka a 10/ Q | J 25 6. i T T T 7 
Q © DSD (T-Matrix) | ° Ka 
> DSD (T-Matrix) fe) DSD | 4 —=Q 
30° > psp (Mie) 7 al | DSD (T-Matrix) 
a 1 ‘i DSD (Mie) : 
L. Q 4 nl 8 
= 25 5 o T 
s 3 | 
o 6- | oy 
= 20} uf 4 @ 4 
2 54 8 
=< g si 
B ist 4 k 
2 15 & 4 
= 104+ , re 
al ; 
A From Wh 
OE ee 0 5 10 15 20 is 
10 15 20 Time (hours UTC) 


Time (hours UTC) 


Time (hours UTC) 
Change in Prediction Error 


2016-04-23 Change in Prediction Error Change in Prediction Error 
T tT t E E 2015-04-27 2016-05-19 
0.5- ° — Change in Prediction Error | | 0.2 t r t t 0.2 tC T : : : 
: No Change ° — Change in Prediction Error * ° — Change in Prediction Error 
0.4; “| 0.15- No Change H 0.15 No Change H 
mo 0.3; t T-matrix More Accurate 4 = 
2 02 0.1; t T-matrix More Accurate 7 0.17 t T-matrix More Accurate 4 
go 0.27 7 
wi 
<c 017 4 4 4 
2 
2 O ere oS ese 2 e ° 
oO e 
o 0.1} 4 
c= 
B -0.2- + a a 
£ 
oO -0.3- 4 Mie More Accurate 4 
“ ; + Mie More Accurate 4 ss 
-0.4- 4 
e 
0.57 1 -0.15 4 4 
r i L r a 
0 5 10 15 20 25 30 0.2 : t t é : 02 : : : : : 
Observed 40 GHz Attenuation (dB) 0 5 10 15 20 25 30 0 5 10 15 20 25 30 


Observed 40 GHz Attenuation (dB) Observed 40 GHz Attenuation (dB) 


Change in Prediction Error 


Change in Prediction Error 
August 2014 - July 2016 


=F : : t 5 tC “T 
e Change in Prediction Error 

BaP 0.5 dB Average i 
t T-matrix More Accurate No Change . 


oO 
S05 
i] 
Ww 
5 0 
© 
i 
a -0.5 
< 
4 
J Mie More Accurate 
-1.5/- 2 
-2 te c c f c cr c cr alk 
5 10 15 20 25 30 35 40 


Observed 40 GHz Attenuation (dB) 


Concluding Remarks 


Conclusions: 


* Prediction of expected attenuation is necessary and may be derived from scaling of existing data at lower frequencies. Use of a 
disdrometer at the ground station site could reduce uncertainties in these predictions of rain attenuation. 


* The Mie and T-Matrix models do not significantly diverge until a droplet diameter of at least 3 mm -- and in the 45° polarization 
case, not until approximately 6 mm or greater, which is rarely observed. 


¢ In many cases, particularly for low rain-rate events, there was found to be little difference between the two models (Fig. 4d), with 
97.2% of observations showing less than 0.1 dB change between the errors of the two predictions, and the change never 
exceeding 1.8 dB. 


¢ At higher rain rates, one model does tend to be preferable but only on a specific case-by-case basis, with a bias toward the Mie 
model being the preferred model more frequently. This may be due to the fact that the disdrometer DSD becomes less 


representative of conditions along the path for larger diameter DSDs, due to its small measurement area. a 
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Thank you! 
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216.433.2546 
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Data Processing Procedure 


Observation Period - Data used in this analysis spans from August 2014 through December 2015 (17 
months) and was processed to isolate applicable rain events. Altogether, 867 hours of rain events were 
isolated across 187 days containing rain. 


Spurious Events - Events were removed from the analysis if the disdrometer was not operational, if the 
receivers were not operational, or if rain occurred along the path but was not measured over the site by 
the disdrometer. 


Attenuation Calibration - For each rain event, the attenuation level before and after the event was 
averaged and subtracted from the measured rain attenuation in an effort to isolate the excess attenuation 
due to rain, though this method is not perfect and other sources of attenuation (e.g. clouds) may still bias 
the result. Instances where this was apparent were removed by inspection (e.g. slow-varying, high 
attenuation consistent with snow or water on the feed/reflector). 


Isolation of Rain Events - Rainy periods were defined as any 1 minute period where the disdrometer 
observed more than 5 precipitating particles (roughly 0.03 mm/hr depending on particle class), and days 
were included in the analysis if they included at least 10 minutes exceeding this level of activity. 


Diameter and Speed Classes 


1 


20.125 0.125 1 2 0.000 0.200 


2 20.250 0.125 2 2 0.200 0.200 
3 20.375 0.125 3 2 0.400 0.200 
4 2 0.500 0.250 4 2 0.600 0.200 
2 0.750 0.250 5 2 0.800 0.200 
6 2 1.000 0.250 6 2 1.000 0.400 
7 2 1.250 0.250 7 2 1.400 0.400 
8 2 1.500 0.250 8 2 1.800 0.400 
9 2 1.750 0.250 9 2 2.200 0.400 
10 2 2.000 0.500 10 2 2.600 0.400 
11 2 2.500 0.500 11 2 3.000 0.400 
12 2 3.000 0.500 12 2 3.400 0.800 
13 2 3.500 0.500 13 2 4.200 0.800 
14 2 4.000 0.500 14 25.000 0.800 
15 2 4.500 0.500 15 25.800 0.800 
16 2 5.000 0.500 16 2 6.600 0.800 
17 25.500 0.500 17 2 7.400 0.800 
18 2 6.000 0.500 18 2 8.200 0.800 
19 2 6.500 0.500 19 2 9.000 1.000 
20 2 7.000 0.500 20 2 10.000 10.000 
21 2 7.500 0.500 . . e 
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Pre-Averaged Error vs. Attenuation 


Instantaneous Scaling Error vs. Attenuation 
August 2014 - December 2015 
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Link Budgets 


19.701 GHz Beacon 


39.402 GHz Beacon 


Parameter User Inputs 


Calculated 


Parameter User Inputs 


Calculated 


Frequency of Operation 19.701 GHz 
Wavelength 
Effective Isotropic Radiated Power (EIRP) 


0.015 m 
19.50 dBW 


Frequency of Operation 39.402 GHz 
Wavelength 
Effective Isotropic Radiated Power (EIRP) 


0.008 
26.50 


Propagation Channel Parameters 


Propagation Channel Parameters 


Transmitter > Receiver Range 38600 km 


Gaseous Absorption Loss 0.5 dB 
Rain Attenuation 0.0 dB 
Pointing Loss 0.0 dB 
Polarization Loss 0.0 dB 
Free Space Loss 


Transmitter — Receiver Range 38600 km 


Gaseous Absorption Loss 0.5 dB 
Rain Attenuation 0.0 dB 
Pointing Loss 0.0 dB 
Polarization Loss 0.0 dB 
Free Space Loss 


Receive Antenna Parameters 

Antenna Diameter 1.2 m 
Illumination Taper Factor 70 deg 
Half Power Beamwidth 
Antenna Efficiency 60 
Antenna Gain 
Noise Temperature Contributions: 

Cosmic Background Noise Temperature 

Atmosphere Physical Temperature 

Antenna Noise Temperature (Clear Sky) 


Antenna Noise Temperature (Rain) 


Receiver Noise Temperature 
System Temperature 


Boltzmann's Constant 

Noise Spectral Density 

Gain over Noise Temperature Ratio (G/T) 
Received Carrier Power (C) 

Carrier to Noise Density (C/NO) 


34.03 
34.03 


634.03 K 
28.02 dBK 
-228.60 dBW/K-Hz 
-200.58 dB 
17.63 dB/K 
-145.41 dBW 
55.17 dBHz 


Receive Antenna Parameters 

Antenna Diameter 0.6 m 
Illumination Taper Factor 70 deg 
Half Power Beamwidth 
Antenna Efficiency 60 
Antenna Gain 
Noise Temperature Contributions: 

Cosmic Background Noise Temperature 

Atmosphere Physical Temperature 

Antenna Noise Temperature (Clear Sky) 

Antenna Noise Temperature (Rain) 

Receiver Noise Temperature 
System Temperature 


Boltzmann's Constant 

Noise Spectral Density 

Gain over Noise Temperature Ratio (G/T) 
Received Carrier Power (C) 

Carrier to Noise Density (C/NO) 


34.03 
34.03 


834.03 K 
29.21 dBK 
-228.60 dBW/K-Hz 
-199.39 dB 
16.44 dB/K 
-144.43 dBW 
54.96 dBHz 
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Hardware Temperature Stability 


LNA Temperature Control (2014-09-03) 
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RF Plate Temperature Control (2014-09-03) 
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IF Enclosure Temperature Control (2014-09-03) 
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System Performance 
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Under normal operating conditions, the 
K-band and Q-band receivers track their 
respective beacon signals independently. 


When attenuation exceeds 30 dB on the 
Q-band channel, the receiver utilizes the 
coherent K-band channel to maintain 
lock on the Q-band (region shown in 
blue). 


Eventually, for deep rain fades, lock can 
no longer be maintained and the noise 
floor of the Q-band receiver is reached. 


Signal lock is immediately regained when 
the signal reappears above the noise 
floor. 
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Measurement Spectral Density 
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Measured Attenuation (2014-06-14) 
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